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INTRODUCTION 

Multiple  mesh  gear-trains  are  very  often  used  In  automotive  transmission 
systems.   The  most  frequently  used  gear-train  is  probably  the  motor  driven 
train  vhich  is  simply  connected  and  is  composed  only  of  spur  gears.   A  schema- 
tical  representation  of  the  n-m.esh  gear-train  system  with  the  notations  to  be 
employed  is  shovm  in  Fig.  1.   The  motor  and  the  load  are  connected  by  the 
pinion-gear  couplings  which,  in  fact,  are  the  vital  part  of  the  study. 

Analysis  of  the  system  reveals  that  for  a  given  motor  torque,  the  accele- 
rations (or  torques)  of  all  shafts  in  the  system  are  entirely  dependent  upon 
the  gear  ratios.   Kence  it  is  desirable  to  determine  the  optimum  gear  ratios 
which  give  the  system  a  maximum  condition  for  a  certain  chosen  criterion.   The 
most  important  criterion  for  the  motor  driven  system  is  the  output  acceleration. 
In  servomechanisms,  the  acceleration  of  the  motor  shaft  is  frequently  desired. 
It  is  to  be  noted  that  for  each  optimization  problem,  the  train  value  or  the 
overall  gear  ratios,  E,  can  either  be  specified  or  unspecified.   Therefore, 
the  chosen  criteria  are:  m.aximization  of 

(i)   the  acceleration  of  the  motor  shaft  with  S  unspecified, 
(ii)   the  acceleration  of  the  m.otor  shaft  with  E  specified, 
(iii)  the  acceleration  of  the  load  shaft  with  E  unspecified,  and 
(iv)   the  acceleration  of  the  load  shaft  with  E  specified. 

The  same  problems  have  been  worked  out  by  the  classical  differential  cal-  • 
cuius  method  for  loss  than  three-mesh  systems.   Burgess  |  1^  presented  the 
solution  of  CASE  (iii)  for  a  one-mesh  system.   For  the  two-mesh  system,  he 
presented  a  sot  of  two  equations  which  are  to  be  solved  simultaneously  to 
obtain  the  optimum  solutions  for  CASES  (iii)  and  (iv).   For  the  n-mesh  system, 
a  recurrence  equation  which  relates  the  optimum  gear  ratios  of  the  successive 
meshes  was  obtained,  and  the  solutions  for  CASES  (iii)  and  (iv)  were  suggested. 
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Ke  also  illustrated  a   method  for  determining  the  optimum  number  of  meshes  by 
means  of  charts.   Mischke  [2]  presented  the  solutions  for  CASES  (i)  and  (iii) 
for  one-mesh  and  two-mesh  systems.  A  comprehensive  inertia  torque  analysis 
was  also  given  by  him. 

The  solutions  of  one-mesh  systems  obtained  by  the  classical  differential 
calculus  method  are  quite  simple.   Whereas  solutions  of  two-mesh,  three-mesh 
systems  are  manageable,  although  a  bit  tedious,  for  a  system  of  more  than 
four-mesh,  the  corresponding  procedures  become  so  involved  as  to  be  practically 
impossible.   Therefore,  the  more  advanced  optimization  techniques  are  needed. 
Here,  a  discrete  version  of  Pontryagin's  maximum  principle  is  used. 

The  maximum  principle  was  first  proposed  in  1956  by  Pontryagin  and  his 
associates  [3~|  for  individual  types  of  tim.e- optimizing  continuous  processes. 
The  first  attempt  to  extend  the  maximum  principle  to  the  optimization  of  stage- 
wise  processes  was  m.ade  by  Rozonoer  [^],  and  various  discrete  versions  of  the 
maximum  principle  were  proposed  by  Chang  [5],  Katz  [6],  Fan  and  Wang  [?!  and 
others. 

The  purpose  of  this  paper  is  to  obtain  the  general  solution  of  n-mesh 
gear-train  systems  for  CASES  (i)  and  (ii)  which  have  not  been  solved,  to  obtain 
a  more  practical  solution  for  CASE  (iii)  and  to  demonstrate  the  advantages  of 
the  discrete  maximum  principle  in  applying  to  the  optimum  designs  of  the  multi- 
stage engineering  systems. 

The  governing  equations  are  developed  first.   The  algorithm  of  the  dis- 
crete maximum  principle  is  introduced,  followed  by  the  solutions  of  n-mesh 
systems  by  this  method.  A  summary  and  two  numerical  examples  are  then  given, 
followed  by  a  conclusion, 

GOVERNING  EQUATIONS  [z] 
The  definition  of  the  gear  ratio  nt  the  n-th  stage  and  its  relation  to 


the  accelerations  of  the  corresponding  shafts  can  be  shown  as  (Fig.  1) 

^      rr 
^„  ITl-l  * 

e"   =  — ^  =  where  m  =  n  (1) 

n 

p       m         * 

The  inertia  forces  analysis  of  a  one-mesh  gear-train  system  can  be 
visualized  by  drawing  a  free  body  diagram  for  each  shaft  (Fig.  2).  For 
shaft  0,  Eulor's  equation  can  be  written  as 

T  -  rV^  =  (I  +  I^)a.  (2) 

M    p  t       M    p  0 

For  shaft  1, 

r-F^  =  (I^  +  r)a,  (3) 

g  t       g    L  1 

where  the  friction  forces,  external  forces  and  inertia  of  the  shafts  are  all 
asstar.ed  to  be  negligible.   Eliminating  F  from  equations  (2)  and  (3)  yields 

1 

\,     =     [iv  +  I^  +-T— (l'  +  ^T^^cc^  (4) 

M     ^  M    P   J.1  ^   g    L  -^ 

^      0 


Since  the  diameters  of  the  gears  in  mesh  are  closely  associated  with  some 

acceptable  contact  ratios,  it  is  often  desirable  that  I  be  replaced  by  some 

g 

equivalent  I  .   For  this  purpose,  it  is  assumed  that  the  gears  in  mesh  are 
made  of  the  same  material  and  have  the  same  thickness,  and  then  the  ratio  of 
the  moment  of  inertias  of  pinion  and  gear  is 

_1       TT  ,  1.4 
^P      -^^ 


1       71         ''     ^ 

g     4  '  g' 


"•'■  The  superscript,  n,  refers  to  the  corresponding  number  of  mesh  counting 
from  the  motor.   The  subscript,  -.,  refers  to  the  corresponding  num.ber  of  shaft 
counting  from  the  motor  shaft. 


shgft_0_./^/ 


shaft    ! 
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Fig.  2.     A     free     body    diagram     of    a     one 
mesh     gear     train     system . 
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or 

S 

■•AT 

(5) 

Substitution  of  th] 

Ls  relation  together  with  equatJ 

Ion  (1) 

into 

equation  (4) 

results 

in 

T 

M 

eq  0  1 

(6) 

where 

"e,^0  '  \* 

1^4 
P 

(e  ) 

(7) 

Equation  (7)  is  regarded  as  the 

equivalent  system 

inertia 

for 

a  one -mesh 

system  referred  to 

the  motor  shaft. 

Expanding  the 

above  analys: 

is. 

it  can  be  proven  that 

the  equivalent  ! 

system 

inertia 

referred  to  the  motor  shaft 

for  a  two-mesh  system 

is 

ec 

I'O   =  C^K-^ 

P 

(e  ) 

(e  ) 

+ 

(e 

P    L-[ 

1  2,2   J 
e  ) 

(8) 

and  for 

each  shaft, 

oq  0  0 

■'m 

oq  1  1 

= 

1   2 
e'(I  )^a, 
cq  0  1 

(9) 

^M 

=   ^4)2^2 

=■ 

1  2, ,2  , 
e  e  (I   )  a„ 
eq  0  2 

In  general,  for  an 

n-mcsh  system,  tho  equivalent  s 

ystem  inertia  referred 

to 

the  notor  shaft  is 


fi"  ^   -  I  ,  iS      P 2.+— —2 E 


(e^S^  +  I^         (e^)S:  -^  I 


p    P  P L         ,,  «x 

,123,2   •^'•'  \  12    n,2         ^^°^ 
(e  G  e  )  (c  G  ...e  ; 


and  for  each  shaft 


T   =  d"  )  a 


M       eq  1  1         eq  0  1 


T„  =  d"  ),a,   =  eV(i^  ).a,  (11) 

eq  2  2  eq  0  <i 


y. 


/T^i  N         1  2     n,_n  . 
T„  =^   (I   )  a    =  e  G  ...  e  (I   )  a 
M       Gq  m  m  eq  U  m 


ALGORITHM  OF  THE  DISCRETE  MAXIMUM  PRINCIPLE  [7,8] 

A  scheinatical  representation  of  a  simple  multistage  system  is  shown  in 

Fig.  3.   The  system  consists  of  N  stage  connGctcd  in  series.   The  state  of 

the  system  stream  denoted  by  an  s-dimensional  vector,  x  =  (x  ,  x  ,  «..,  x  ), 

1   <i        s 

is  transformed  at  each  stage  according  to  an  r-dimensional  decision  vector, 
0  =  (0  ,  0„,  •••♦  0  ),  which  rcpresonts  the  decision  made  at  that  stage. 
The  transformation  of  the  system  stream  at  the  n-th  stage  is  described  by  n 
set  of  performance  equations. 
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Fig.  3.    Multistage       decision        system.    . 
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n-1  ^n 

(12) 

X. 

1 

=  a_ 

■-   1,2,-. -,3; 

n  =  1, 

,2,.. . ,N 

or  in  vector  form 

n 

X 

=  T^(x"-^;  e")   , 

n  =  1,2,. 

..,N 

(13) 

0 

X 

=  a  . 

A  typical 

.  optimization  problem 

associated  with 

such  i 

a  system 

is  to  find 

a  sequoncG  of 

e". 

n  =  1,2,...,N,  subject  to  constraints 

f1 

[e^  ej,  ...,  9^] 

n 
<  0  , 

i 

=  1,2 
=  1,2 

,...,N 
,...,r 

(14) 

which  nakes  a 

function  of  the  state 

variable  of  the 

final 

stage 

S 

=  i;   c.x.  , 
i=i    '  ' 

c  =  constant 

(15) 

an  extreTTUKi  when  t 

:he  initial  condit: 

0 
Ion  X  =  a  is  given. 

The  function  S  which 

is  to  be  naxii; 

:ized 

'.    (or  niniTnized) ,  ; 

is  the  objective 

funct 

ion  of  the  system. 

The  procedure 

■  for  solving  such 

an  optimization 

probl 

em  by  a  c 

liscrete 

version  of  the  maximuin  principle  is 

to  introduce  an 

s-dimensional 

adjoint 

n 
vector  z  and 

a  Hamiltonian  function  H  which  satis' 

fy  the 

following  relations: 

h'^ 

=  (z'')V  =  Z  zV(x"-^;  e")  , 

i=l  i  i 

n 

=  1,2,... 

•,N     (16) 

n- 
z 
1 

■1     S  k" 

>  n-1 
ox.   , 

i  =  1,2,...,, 
n  =  1,2,... , 

s; 

(17) 

and 

N 

=   C^  ,          i  : 

=  1,2,. ..,N  . 

(18) 

10 

If  the  optimal  decision  vector  function,  9  ,  which  makes  the  objective 
function  S  an  extremum  (maximum  or  minimum),  is  interior  to  the  set  of  admis- 
sible decision,  e""",  a  necessary  condition  for  S  to  be  a  (local)  extremum  with 
respect  to  0   is 

li^  =  0  ,       n  =  1,2,. ..,N  (19) 

If  ?^  is  at  the  boundary  of  the  set,  it  can  be  determined  from  the  condition 
that  H'  is  (locally)  extremum. 

SOLUTION  BY  THE  DISCRETE  MAXIMUM  PRINCIPLE 

The  optimum  solutions  for  the  four  chosen  criteria  can  be  obtained  by 
formulating  the  equations  derived  previously  to  fit  into  the  algorithm. 

(i)  Maximization  of  the  acceleration  of  the  motor  shaft,  a^,  with  the  train 

value,  E,  unspecified. 

An  N-mesh  gear-train  system  is  represented  schematically  in  Fig.  4,  where 
a  stage  represents  a  mesh  and  the  n-th  stage  represents  the  n-th  mesh  counting 
from  the  motor.   The  decision  variable  0  and  the  state  variables  x^  and  x^ 
are  defined  as  follows: 

n  n    n 

9  =  the  gear  ratio  at  the  n-th  stage,  hence  9   =  e   , 

X  =  the  accumulated  train  value  up  to  and  including  the  n-th  stage 

of  r;-stage  system, 
X  ^  the  equivalent  system  inertia  referred  to  the  motor  shaft  up 

to  and  including  the  n-th  stage. 

Thus  for  the  first  state  variable,  the  performance  equation  is 

.";  =  x;;-v,      n.i,2,...,N 

(20) 
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For  tho  second  varia'rle;  for  n=l,  from  ecu?.*: .'.on  (7),  x^  ccn  be  ^Titter. 

as 

x"  =  -  +  '  +  ~ 

■2            -M   -P      ^^1,2 

(21) 

if  X  is  defined  as 

x°  =  I  -J 

2      M    L 

this  equation  is  substituted  in  -qv.nt'.on  ("1)  r.nd  results  in 

x^  =  X  + ■■—r-T 

(22) 

2 
For  n=2,  from  equation  (8),  x  can   bo  vrltten  as 

^2     =  K   "^  ^        1  -   ~  "^     12  2 

2        >.      ?           (gX)'-              (e-s^)"^ 

(23) 

Substitution  of  equation  (22)  into  (2:.)  results  in 

.^r.^f  ..  (r^   .  T  ^rs2)2  ^i 

(24) 

and  it  can  be  shoTv'n  that  "e  cr.n  ^■T:r-ite  x^ ,  in  general,  as 

n      n-1    -T^-  '   +  '^^•^   -■  -T^-'^y  '   ^  \ 
X   =  x^   4- -■ 

(25) 

The  perfomance  equctic-'.s  fcr  the  probler,  are  si;:rjnariz;ed  as  follows. 

jt      r  - 1  n 
>--,   =  •<■  9   ,       n  .-^  :,2,...,N  , 

X            i. 

(26) 

0 

X   =  1 

1 

13 

n     n-1      y              ,  o     m 

^"2  -    \       ^'     n-1  n2        ""       ''"'   '"^  ' 
(x,  9  )  ♦ 

(27) 

x°  =  I  +  I 
2      ML 

where 

y"    .    i^oVh-   (I^-  I  )(e")^   +  1                (28) 

p         p    L         L 

subject  to  constraints 

e^  >  1  ,      n  =  1,2,. ..,N                          (29) 

The  objective 

function  is  to  riinimize 

2     K      N 
3  =  1  C.X.   =  X   ,                                 (30) 

i=l   ^  ^      ^ 

• 

c,  =  0,       c^  =  1 

which  is  the 

total  equivalent  inertia  referred  to  the  motor  shaft  of  the  N 

stage  system. 

The  Hamiltonian  function  and  the  adjoint  vectors  are  as  follows. 

H^ 

np  n-l.n-i    nr  n-1       v"   i           «        .  x 
=  z^[x^  e  ]  +Z2L-2   ■"       n-ln2^'     n^l,2,...,N   (31) 

(x^   0  ) 

n-1      Sh"       n„n       2  (y^) 

i.       N  n-i      1       n-l.j>,^n.^ 
^x^           (x^  )   (e   )  . 

z'^   -  c^  =  0                                       (33) 

n-1      ^  H        n 
^2    ~  V  n  1  "  ^^2  '       "  -  l'2,-*-»^'             (34) 
^^2 

N 
^2  =  ^2  "  ^                                       (35) 

14 


Fron  equations  (34)  and  (35),  it  concludes  that 


z^  =  1,       n  =  1,2,. -.^N  (36) 


Hence   the   Haniltonian  function,  equation   (31),    is  reduced  to 


,n  n  n-l.n    ,     n-1    ,  y  i    o  >•  /o-7> 

H       =     zx,      9     +  x^       + :; ^  n  =  1,2,...,^  (37) 

11  2  n-i.  Ti  I 

(x,      9   )      » 


n  n 

Partially  differentiating  H     with  respect  to  9     gives 


or 


^9"       '■'■       (x;^-b^9^)^ 


^           .  2[i;(9")^  -    I    ] 
z       =     ■ —      n  =  l,2,...,N  (38) 


Substituting  equation  (38)  into  both  sides  of  equation  (32)  gives  the  recurrence 
relation 

lj"^e"-^^  =  Ip[2(e")'-M],       n  =  2,3,...,N     (39) 

which  relates  the  optiinu.-n  gear  ratios  of  two  successive  gear  ir.eshes.   Reaem- 

N 
bering  that  z     =  0,  equation  (38)  can  be  written  for  r.  =  N,  as 

1^9^')'  -  I,  =  0 

P        L 

or 

N        I.  k 
^        =      [41  (40) 

I" 
P 

Substituting  9   into  the  recurrence  relation,  equation  (39),  9  "  can  be 
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obtained.   Likowisc  9^^,  n  =  N-2,  N-3,  •••,  1  can  be  determined  and  the  set  of 


9   is  the  optimum  solution. 


For  instance,  for  a  one-mesh  gear-train  system  (N=l),  from  equation  (40), 


This  agrees  with  the  result  obtained  by  Mischke  [2],   For  a  two-mesh  gear- 
train  system  (N=2), 

^   -  L  ^  J 
I 
P 


e 


2   o  /  2 
1       I  +  2  /  I  I^   ^ 
'■  r    P         V  p  L 


I 
P 


P  iv-] 

1       -' 


The  optimum  train  value  can  be  calculated  from  equation  (26).   The  minimum 

N 
equivalent  system  inertia  referred  to  the  motor  shaft,  x^ ,  can  be  obtained  by 

equation  (27),  and  the  maximum  acceleration  of  the  motor  shaft,a-,,  is  to  be 

determined  from 

(ii)  Maximization  of  the  acceleration  of  the  motor  shaft,  a,.,  with  the  train 
value,  E,  specified. 
Since  the  train  value,  E,  is  specified,  the  problem  is  usually  classified 
as  a  fixed  end  problem  and  the  aim  *.s  to  apportion  E  among  the  N  stages  so 

that  a-|  is  a  maximum.   The  definitions  of  the  decision  variable  9  and  the 

.  ,  ,    ri      n 
state  variaoles  x  and  x^  remain  the  same  as  in  case  (i),  while  the  performance 

equations,  the  constraints  and  the  objective  function  arc  also  remained  un- 
changed except  that 
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is  to  be  added  to  equation  (26),   Hence  instead  of  equation  (33), 

z^'  ^  c   =0  («) 

By  following  the  sair.e  procedures,  it  ends  up  with  the  same  recurrence  relation 
as  for  CASE  (i).   This  relation  is  repeated  as  follows 

jn-l(gn-1^4  ^  i^r2(9")^  +1]  ,      n  =  2,3,...,N     (44) 
P  P^ 

But  for  this  case,  an  equation  equivalent  to  equation  (40)  of  CASE  (i)  is  not 
available.   Nevertheless,  the  set  of  optimum  0  can  be  uniquely  determined  by 
equations  (44)  and  (42).   The  solution  can  bo  obtained  numerically  and  will 
be  illustrated  in  the  next  section.   Once  the  set  of  9   is  known,  it  is  pos- 
sible to  calculate  x_,  hence  a„  as  mentioned  in  CASE  (i). 

(iii)  Maximization  of  the  acceleratior  of  the  load  shaft,  a  ,  with  the  train 

m 

value,  E,  unspecified. 

The  nature  of  the  problem  remains  the  same  as  in  CASS  (i).   Instead  of 

maximizing  the  acceleration,  the  equivalent  system  inertia  referred  to  the 

load  shaft,  (I   )  ,  shall  be  minimized.   The  decision  and  state  variables  are 
eq  m 

defined  as  follows 

6  =  the  gear  ratio  at  the  n-th  stage,  hence  9  =  e  , 

X  =  the  accumulated  train  value  up  to  and  including  the  n-th  stage 


of  N  stage  system, 
'2 


n 
X  =  the  equivalent  system  inertia  referred  to  the  load  shaft  up  to 


and  including  the  n-th  stage. 
Thus  for  the  first  state  variable,  the  performance  equation  is 
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n 

X   e   ,       n  =  1,2,...  ,N 

(45) 

0 

1 

For  the  second  state 

1 
variable:  for  a  one-stage  system,  n=l,  x"  can  be 

written 

as 

1 
X2  = 

1    ^>'>'-^T 

(46) 

again  defining 

0 

^2  = 

then  substituting  this  relation  into  equation  (45)  results  in 

1 
^2  = 

O^X    ^>'>'-^<Ip-  V^^'^'-^'^L 

(47) 

2 
For  n=2,  from  equations  (8)  and  (9),  x  can  be  written  as 

2 
^2  = 

1  ^y^'-^i  ^y^'-h 

(4S) 

Substituting  equation 

(47)  into  (48)  yields 

2 

2   2  4     2        2  "> 

(49) 

'^    ^ 

and,  in  general,  x  c; 

in   be  written  as 

n 

"2 

.n-l.n  .  I>V  .  (I^  V(eV  .  I, 

(50) 

^                    e'o'...e" 

n  -  1,2,. ..,x 
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The  performance  equations  are  suiranarized  as  follows 

n      n-l^n           ,  ^ 
x,^  =  x^  9   ,       n  =  1,2,. ..,N 

(51) 

0 
X    =   1 

1 

1        n 
n      n-l^n     y              ,  ^     ^. 

^"2   -  "2   ^   "^  n-l  n         "   ^'^'    '^ 

(52) 
0 
X   =1+1 
2      M    L 

subjects  to  constraints 

0^  ^  1  ,       n  =  1,2,-..,N                         (53) 

The  objective  function  is  to  iriininize  the  total  equivalent  system  inertia 

referred  to  the  load  shaft 

t             ^               ^ 
S  =  Z  ex.  =  X 

i=l   ^  ^      2 

(54) 

c^  =  0  ,       c^  =  1 

The  Hamiltonian  function  and  the  adjoint  variables  are 

h"  -  z"(v''-^e'')  +r"r^^-io"  ^,  y"  -   „   1  -  ...  V  (55) 

H      z^l..^  ""^^^■I'-'^l      °   ■   n-l  n-  ,  n  -  l»-,---,N  (5j) 

X   0 
1 

n                    ri 
n-l       c^  H        n^n    Hp     y     , 

^    ^  .  n-l  -     \^      -\^-.       .-I2n^*              (5^> 
^^1                 ^   ^^     n  =  l,2,...,. 

\      -  c^  =  0                                       (57) 
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n-1      h   H       n^n            ,  _     ., 

z_    = =  z  0   ,       n  =  l,2,..',i\ 

2       ^  n-1      ^ 
o  x^ 

(58) 

^2  -  c^  =  1 

(59) 

Equations  (58)  and  (59)  gives 

• 

n      n+l„n+l            ,  _     .,  , 
Z2  =  Z2  ®    '       "  "^  1,2,«",N-1 

(60) 

Taking  the  partial  derivative  of  equation  (55)  with  respect  to  9   yields 

or 

i 

np^  n,^n-4   ,  n-1  n-1   ^n   ^  s^rP-s."^        -r  "> 

(61) 

'1  "                  <)' 

n  =  1,2,...,N 

Substitution  of  equation  (61)  into  both  sides  of  (56)  and  making  use  of 

equation  (60)  result  in  the  following  recurrence  relation 

in-l(gn-1^4.^  ^n^2(9V  +1I,       n  =  2,3,...,N 
P              P 

(62) 

which  is  identical  to  equations  (39)  and  (44),  the  recurrence  relations 

J 

1 

obtained  for  CASES  (i)  and  (ii).   For  n=N,  equation  (60)  becomes 

1 

P         1   2      p    L         L 

(63) 

Thus  for  a  one-sta^e  systein,  N=:l,  equation  (63)  becomes 

■ 

3lJ;(0^^  +  (T  4-  i^o')^  -1^0 
P          M    p          L 

(64) 

i 
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The  solution  of  this  equation  gives 


-  (I  +  i^)  +  /(I  +  i^)^  +  121^1     h 
e^  =  [ ^- E V  ^   P EJi]        (65) 

which  agrees  with  the  solutions  obtained  by  Burgess  [ij  and  that  by  Mischke  [2], 

n 
However,  for  N  equals  two  or  more,  the  set  of  optintm  9   is  to  be  determined 

by  equations  (62)  and  (63).   The  solution  can  be  obtained  numerically  and  will 

be  illustrated  in  the  next  section.   Once  the  set  of  9   is  known,  it  is  pos- 

N 
sible  to  calculate  x_,  and  hence  a  by 

i.  ra 

T,  =  x!^a  (66) 

M      2  m 

(iv)  Maximization  of  the  acceleration  cf  the  load  shaft,  a.  ,  with  the  train 

m 

value,  S,  specified. 
As  in  CASE  (ii),  the  problem  is  classified  as  a  fixed  end  problem  and  the 
aim  is  to  apportion  E  among  the  N  stages  so  that  a   is  a  maximum.   The  defi- 
nitions of  the  decision  variable  9  and  the  state  variables  x  and  x'  remain 
the  same  as  in  CASE  (iii),  while  the  performance  equations,  the  constraints 
and  the  objective  function  are  also  unchanged  except  that 

N 
x^  =  E  (67) 

is  to  be  added  to  equation  (51).   Hence  instead  of  equation  (57), 

By  following  the  same  procedures,  it  ends  up  with  the  following  rocurrcnce 
rclat Ion 

l';-^0"-^)'^  .   l"[2(0")'  4-1]  ,       n=.2,3,...,N    (69) 
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which  is  identical  to  equations  (39),  (44)  and  (62),  the  recurrence  relation 
obtained  for  CASES  (i),  (ii)  and  (iii).   The  opt-'Tnutn  solution  9  can  be  de- 
termined by  equations  (67)  and  (69).   It  is  to  be  noted  that,  in  spite  of  the 
differences  betveen  the  performance  equations,  the  tvo  equations  which  de- 
termine the  optimum  solution  for  CASE  (iv)  are  identical  to  those  of  CASE  (ii). 
Therefore  CASE  (ii)  and  CASE  (iv)  are  called  the  dual  problem,  i.e.,  to  maxi- 
mize CASE  (ii)  is  equivalent  to  maximize  CASE  (iv). 

SUMMARY  AND  NUMERICAL  EXAMPLES 

In  spite  of  the  differences  in  the  performance  equations  and  the  criteria 
chosen  to  be  m.aximized,  it  has  been  found  that  the  recurrence  relations  which 
relate  the  gear  ratio  of  the  successive  gear  mesh  are  identical  for  the  four 
cases  considered.   This  relation  is  repeated  below. 

jn-l^^n-1^4  ^  i^p^oV+l]  ,       n  =  2,3,...,N    (70) 

which  can  be  re'snritten  as 

n^        n  2 


n-x 


=     r_E 


io[2(e')    +1]^^. 


-1 ]     ♦  n  =  2,3,. ..,N  (70a) 

n-1 


^P 


or 


e       =      [ :; -^    ,  n  =  2,3,.",N  (70b) 

This  recurrence  relation  consists  of  N  unknowns  but  with  only  N-1  equations. 
An  additional  equation  obtained  for  each  of  CASES  (i)  and  (iii)  and  given  for 
each  of  CAS!^S  (ii)  and  (iv)  is  rewritten  here: 
CASE  (i) 

^   "   I--7tJ  (40) 

r 
p 
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CASE  (ii) 

, 

N 
X   =  E     (given)                                 (^2)      ; 

CASE  (iii) 

■1 

31  (9  )   +  (x,   X.   +  I  -  I  )(9  )   -  I   =  0           (63) 

p         1   2      p    L         L                         } 

CASE  (iv) 

N 
x^   =:  E     (given)                                  (67) 

Solution  of  CASE  (i)  is  quite  straightf onward,  while  the  solutions  of 

CASES  (ii). 

(iii)  and  (iv)  are  to  be  obtained  by  numerical  iteration  methods. 

A  specific  one,  the  Falsi  method,  has  been  applied  to  obtain  the  solution  as 

described  in 

the  steps  below. 

CASE  (ii) 

Step  1. 

Assume  two  trial  values  for  9  ,  i.e.,  9  (1),  9  (2).                 i 

Step  2. 

Calculate  9  ,  n  =  2,3,*",N  by  means  of  equation  (70b).             j 

Step  3. 

Calculate  x  ,  n  =  1,2,«»»,N  by  means  of  equation  (26). 

Step  4. 

Calculate 

N                                                        '■ 

G   =:   x^  -  E 

Step  5. 

Check  the  inequality                                            ' 

G     -   ^     ^    0                                                                                                                            1 

whore  (S  is  the  maximum  allowable  error.   If  it  is  satisfied. 

the  sot  of  e   is  the  optimum  solution.   Go  to  Step  10.               j 

Step  6. 

Otherwise  proceed  to  the  next  step. 

Calculate  a  now  trial  value  9  (3)  from  the  following  equation. 

• 

0^3)   .  ^'(^^^(^^     -    0^1)  G(2)  ■ 

r.(i)  -  G(2)                      ^'^^ 
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Step  7.    Repeat  Stop  2  through  Step  4  to  obtain  G(3). 
Step  8.    Check  the  inequality 

|g(3)|  -6^0 

If  it  is  satisfied,  the  set  of  G   is  the  optimum  solution. 

Go  to  Step  10.   Othem-rise,  proceed  to  the  next  step. 
Step  9.    Replace  9  (2)  and  0~(3)  by  9  (1)  and  9^(2),  and  G(2)  and  G(3) 

by  G(l)  and  G(2),  respectively.   Repeat  Step  6  through  Step  8, 
Step  10.   Calculate  x  ,  n  =  1,2,...,N  by  means  of  equation  (27). 
Step  11.   Calculate  a»  by  means  of  equation  (41). 

CASE  (iii)   The  procedures  given  in  CASE  (ii)  can  be  applied  to  this  case  by 
making  the  following  modifications. 

(1)  Replace  Step  4  by  Steps  4a  and  4b  shown  below. 

Step  4a.   Calculate  x  ,  n  -   1,2,'««,N  by  means  of  equation  (52). 
Step  4b.   Calculate 

G    =   3  AeV  +  (xf-^x^'-^  H- i^^  -  T)(eV-  I, 

p  1    ^       p     L  L 

(2)  Replace  Step  10  by:      Calculate  a     by  means   of  equation   (66). 

m 

(3)  Skip  Step  11. 

CASE  (iv)   Since  CASES  (iv)  and  (ii)  are  dual  problems,  the  procedures  of  solu- 
tion are  identical  to  that  of  CASE  (ii). 

The  determination  of  the  optimum  number  of  meshes  for  a  certain  chosen 
criterion  can  be  done  by  tabulating  and  comparing  the  optimum  solutions  obtained 
for  n  =  1,2, • •• ,n. 

The  following  numerical  examples  are  solved  as  an  illustration. 

Example  1.   The  optim.um  gear  ratio  of  each  mesh  for  m.aximizing  the  four  chosen 
criteria  for  a  throe-mosh  and  n  five-m.esh  gear-train  system  is  obtained.   The 
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given  parameters  for  these  systems  are  shovm  in  Tables  1  and  2. 

Example  2.   The  parameters,  I^^,,  I,^,  T^,,,  I^  and  I^  of  an  existing  two-mesh 
gear-train  system  are  given  in  Table  3.  Additional  pinions  with  I^  =  0.3 
are  available.   The  optimum  additional  n'jmber  of  gear  m,eshes  to  be  installed 
in  the  system  for  maximizing  the  acceleration  of  the  load  shaft  with  the  train 
value  unspecified  is  obtained. 

The  solutions  of  Examples  1  and  2  have  been  obtained  by  the  use  of  a  com- 
puter and  are  tabulated  in  Tables  1,  2  and  3.   The  results  of  simulations  indi- 
cate that  the  solutions  obtained  are  actually  the  optimum  ones. 

Referring  to  CASE  (i)  and  CASE  (ii)  in  Tables  1  and  2,  for  the  maximization 
of  the  acceleration  of  motor  shaft,  a^,  with  either  the  train  value,  E,  speci- 
fied or  unspecified,  the  optimum  gear  ratios,  9  ,  are  allocated  in  such  a  way 
that  the  first  N-1  gear  ratios  are  close  to  one  another,  while  the  gear  ratio 
at  the  last  stage,  G",  is  m.uch  larger  which,  in  fact,  is  the  decisive  factor 
for  the  system  train  value,  E.   The  maximum  acceleration,  a^^,  is  almost  inde- 
pendent of  the  train-value,  E,  in  other  words,  no  matter  what  train  values  are 
assigned  to  the  system,  the  maximum  motor  shaft  acceleration  remains  almost 
the  same.   For  CASE  (iii),  the  maximization  of  the  acceleration  of  load  shaft, 

a  ,  with  the  train  value,  E,  unspecified,  the  gear  ratios,  9  ,  are  all  close 

m 

to  one  another,  and  there  is  no  decisive  factor  for  the  train  value.   While  for 

CASE  (iv),  with  the  train  value,  E,  specified,  the  gear  ratio  of  the  last  mesh, 

K 
9  ,  is  the  decisive  factor  of  the  train  value  as  appeared  in  CASES  (i)  and  (ii). 

Reforrirg  to  Table  3,  the  results  show  that  a  four -mesh  system  gives  the 

largest  maximum  load  shaft  acceleration,  a  •   Thus,  two  more  meshes  are  to  be 

m 

installed  in  the  system.  However,  it  must  be  pointed  out  that  this  optimal 
policy  is  for  maximization  of  the  load  shaft  acceleration  as  the  criterion. 
It  might  not  bo  the  optimal  policy  for  the  system  if  cost  factor  is  considered. 
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The  examples  have  been  worked  out  merely  to  illustrate  the  procedures  of 
solutions.   Therefore,  no  dimension  has  been  assigned,  and  the  given  parameters 
have  been  arbitrarily  chosen.   However,  the  inertia  of  the  pinion  chosen  must 
always  satisfy  the  following  relation  which  is  derived  from  equation  (70b). 

jn-l(Qn.lj4  ^  jn 
P  P 

CONCLUSION 

By  the  use  of  the  discrete  maximum  principle  with  the  aid  of  modern  com- 
puters, the  iJractical  optimum  solutions  of  n-mesh  gear-train  systems  for  the 
four  chosen  criteria  and  the  optimum  number  of  meshes  can  be  determined.   The 
study  reveals  that  the  recurrence  relations  which  relate  the  optimum  gear-ratios 
of  the  gear-train  systems  are  identical  for  all  four  cases.   For  optimization 
of  the  acceleration  of  any  intermediate  shaft  between  the  motor  and  the  load 
shafts  when  the  train-value  is  specified,  the  recurrence  relations  are  also 
identical  to  the  one  obtained,  whereas  for  the  case  that  the  train  value  is 
unspecified,  the  recurrence  relation  remains  to  be  investigated.   The  feature 
of  the  discrete  maximum  principle  is  that  the  n-stage  system  as  a  whole  is 
treated,  and  hence  the  solutions  obtained  are  valid  for  the  n-stage  system  in 
general,  and  the  computational  procedures  are  simpler  than  the  conventional 
methods.   It  shows  that  the  discrete  maximum  principle  is  a  powerful  tool  in 
dealing  with  the  optimal  designs  of  multistage  systems.   It  can  be  applied  to 
obtain  the  optimum  solutions  for  many  other  design  problems. 


29 


REFERENCHIS 

1.  Burgess,  Jr.,  E.  G. ,  "Minimization  of  Goar-Train  Inertia,"  Trans.  ASME, 
76,  493-496  (1954) 

2.  Kischke,  C.  R. ,  "lements  of  Mechanical  Analysis,  Addison-Wesley  Publishing 
Coinpany,  Inc.,  Reading,  Mass.,  (1963). 

3.  Pontryagin,  L.  S.,  V.  G,  Boltyanskii,  R.  V.  Gamkreldze,  and  E.  F.  Mischenko, 
The  Mathematical  Theory  of  Optimal  Process,  (English  translation  by 

K.  N.  Trirogoff),  Interscience,  New  York,  (1962). 

4.  Rozonoer,  L.  I.,  "The  Maxinum  Principle  of  L.  S.  Pontryagin  in  Optimal- 
Systen  Theory,"  Automat,  Telmech.,  Moscow,  2Ci,  1320,  1441,  1561,  (1950). 

5.  Chang,  S.  S.  L.,  "Digitized  Maximum  Principle,"  Proc.  IP^,  4S,  2030-2031 
(1960). 

6.  Katz,  S.,  "A  Discrete  Version  of  Pontryagin's  Maximum  Principle,"  J. 
Electron.  Control,  13,  179  (1962). 

7.  Fan,  L.  T. ,  and  C.  S.  Wang,  The  Discrete  Maximum  Principle  —  A  Study  of 
Multistare  Systems  Optimization,  Wiley,  New  York, (1964). 

8.  Fan,  L.  T.,  C.  L.  IVang,  and  C.  S.  Wang,  "Optimization  of  Step  Rockets  by 
the  Discrete  Maximum.  Principle,"  Trans,  of  ASME,  J.  of  Engineering  for 
Industry,  ST^,  Ser.  B,  419-424  (1965). 


30 


^X)MENCL^TURE 

e  =  gear  ratio  or  step  dovn  ratio  oi:  nuh  stage 

E  =  overall  gear  ratio  or  train  value 

F  =  transmission  force 

H  =  Hamiltonian  function,  defined  by  equation  (16) 

d"  )   =  equivalent  system  inertia  of  n-mesh  system  referred  to  mth  shaft 
eq  m 

I   =  inertia  of  gear  of  nth  stage 
g 

T   =  inertia  of  load 

I   =  inertia  of  motor 

M 

I   =  inertia  of  Pinion  of  nth  stage 
P 

N  =  number  of  stages  of  a  raultim.esh  gear- train 

n  =  the  nth  stage  of  a  multimesh  gear-train  defined  in  Fig.  3 

r  =  nunber  of  decision  variables 

r  =  radius  of  gear  of  nth  stage 

or 
O 

r   =  radius  of  pinion  of  nth  stage 
P 

S  =  objective  function 

s  =  number  of  state  variables 

T  =  transformation  poerator,  defined  by  equation  (12) 

T   =  motor  torque 

K 

X  =     state  vector 

y  -     function  defined  by  equation  (2S) 

z  =  covariant  vector  of  x  introduced  in  maximum  principle 

a   =  acceleration  of  the  mth  shaft 
m 

^    -     constraint  function 
9   =  control  variable  of  nth  stage 
Superscript  n  =  stage  nuirbor 
Subr.cript  m  -  shaft  number 
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The  general  algorithm  of  the  discrete  maxinuni  principle  is  stated  briefly. 
It  is  applied  to  obtain  the  optimum  solutions  of  design  problems  involving 
the  n-T^esh  gear-train  system  for  four  chosen  criteria  which  are:  the  maximiza- 
tion of  i)  the  acceleration  of  the  motor  shaft  with  the  train  value  unspecified, 
ii)  the  acceleration  of  the  motor  shaft  with  the  train  value  specified,  iii) 
the  acceleration  of  the  load  shaft  with  the  train  value  unspecified  and  iv) 
the  acceleration  of  the  load  shaft  with  the  train  value  specified. 

The  purpose  of  this  report  is  to  obtain  the  general  solution  of  n-mesh 
gear-train  systems  for  the  chosen  criteria  (i)  and  (ii)  which  have  not  been 
solved,  to  obtain  a  more  practical  solution  for  the  chosen  criterion  (iii)  and 
to  demonstrate  the  advantages  of  the  discrete  maximum  principle  in  applying 
to  the  optim.um  designs  of  the  multistage  engineering  systems. 

The  study  reveals  that  by  the  use  of  the  discrete  maximim  principle  with 
the  aid  of  modern  computers,  the  practical  solutions  of  optimization  problems 
involving  design  of  n-mcsh  gear-train  system.s  for  the  four  chosen  criteria  and 
the  optimum  number  of  m.eshes  can  be  determined.   It  also  reveals  that  the 
recurrence  relations  which  relate  the  optimum:  gear-ratios  of  the  gear-train 
systems  are  identical  for  all  four  cases.   The  discrete  maxim'^m  principle  is 
shown  to  be  a  powerful  tool  in  dealing  with  the  optim.al  designs  of  multistage 
systems.   It  can  be  applied  to  obtain  the  optim.um  solutions  for  many  other 
multistage  design  problem.s. 


